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Abstract. We show that the 3D wedge filling transition in the presence of short- 
ranged interactions can be first-order or second order depending on the strength 
of the line tension associated with to the wedge bottom. This fact implies the 
existence of a tricritical point characterized by a short-distance expansion which 
differs from the usual continuous filling transition. Our analysis is based on an 
effective one-dimensional model for the 3D wedge filling which arises from the 
identification of the breather modes as the only relevant interfacial fluctuations. 
From such analysis we find a correspondence between continuous 3D filling at 
bulk coexistence and 2D wetting transitions with random-bond disorder. 
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Figure 1. Schematic illustration of a typical interfacial configuration and relevant 
lengthscales for a fluid adsorption in a 3D wedge. 

Fluid adsorption in micropatterned and sculpted geometries has become the 
subject of intense study over the last decade. Highly impressive technological advances 
which allow the tailoring of micro-patterned and structured solid surfaces on the 
nanometer to micrometer scale PP are a landmark in the development of the emerging 
microfluidic industry |5j which aims at minituarizing chemical synthesis plants or 
biological analysis equipment in much the same way the silicon chip brought about 
the electronics revolution. However, the theoretical understanding of this phenomenon 
is far from being complete. Recent studies of filling transitions for fluids in 3D wedges 
show that interfacial fluctuations are greatly enhanced compared with wetting at flat 
substrates 00]. The control of such enhanced interfacial fluctuations is crucial for 
the effectiveness of the microfluidic devices. Fortunately, there are simple theoretical 
approaches which take into account these effects. For example effective Hamiltonian 
predictions for the critical exponents at continuous (critical) wedge filling with short- 
ranged forces have been confirmed in large scale Ising model simulation studies p]. 
Similar experimental verification of the predicted geometry-dominated adsorption 
isotherms at complete wedge filling raise hopes that the filling transition itself and 
related fluctuation effects will be observable in the laboratory. We further develop 
the theory of wedge filling in this paper, focussing on the emergence of a new type of 
continuous filling: tricritical filling. 

First we briefly review the fluctuation theory of 3D wedge filling. Consider the 
interface between a bulk vapour at temperature T and saturation pressure with a 
3D wedge characterised by a tilt angle a. Macroscopic arguments dictate that the 
wedge is partially filled by liquid if the contact angle 9 > a and completely filled if 
9 < a 0. The filling transition refers to the change from microscopic to macroscopic 
liquid adsorption as T — > Tf, at which 9(Tf) = a, and may be first-order or continuous 
(critical filling). Both of these transitions can be viewed as the unbinding of the liquid- 
vapour interface from the wedge bottom. Characteristic length scales are the mean 
interfacial height above the wedge bottom lw, the roughness £j_ and the longitudinal 
correlation length £ y , measuring fluctuations along the wedge (see Figure 0). The 
relevant scaling fields at critical filling are 9 — a and the bulk ordering field h (which 
is proportional to the pressure difference with the saturation value). At coexistence 
(h = 0) we define critical exponents by lw ~ (0 — a)~^ w and £ y ~ (6 — a)~ Uy . The 
roughness can be related to £ y by the scaling relationship £j_ ~ Q w , where C,w is the 
wedge wandering exponent. For short-ranged forces, £w = 1/3. 

For shallow wedges, i.e. a <C 1, the free energy of an interfacial configuration can 
be modelled by an effective Hamiltonian based on the capillary wave model for wetting 
of planar substrates However, an analysis of this model |3] shows that the liquid- 
vapour interface across the wedge is aproximately flat and soft-mode fluctuations arise 
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from local translations in the height of the filled region along the wedge. These breather 
modes are the only relevant fluctuations in the continuous filling phenomena, and can 
be taken into account by the following effective Hamiltonian [3] 



where lo(y) is the local height of the interface at position y along the wedge bottom 
and S is the liquid-vapor surface tension. Note that the effective bending term 
resisting fluctuations along the wedge is proportional to the local interfacial height. 
The effective binding potential Vw(lo) f° r h ^ lin where 1^ is the mean wetting film 
thickness for a planar substrate, is given (up to irrelevant additive constants) by |3J: 



where W(l) is the binding potential between the gas-liquid interface and a planar 
substrate. Note that the mean field result for lw is recovered by minimizing Vw{lo) 
for Iq > l„, which is an important check on the self-consistency of the method PJ- In 
addition to a hard wall repulsion for Iq < 0, the potential Vw(h) contains a short- 
ranged attraction which may be modified by micropatterning a stripe along the wedge 
bottom, so as to weaken the local wall-fluid substrate and therefore strengthen the 
interfacial binding. This observation will be crucial for the existence of tricritical 
filling, since with this decoration it may be possible to bind the interface to the wedge 
bottom at the filling boundary 8 — a and h = 0. For latter convenience, hereafter we 
will set h = in our discussion of continuous filling. 

The quasi-one-dimensional nature of the effective Hamiltonian Equation l|T|l 
allows us to use the transfer-matrix formalism. In the continuum limit the partition 
function is defined as a path integral J2j (setting fc^T = 1 for convenience) 



where Y is the wedge length and l a and lb are the endpoint heights. The position- 
dependent stiffness introduces some ambiguity in the definition of the path integral. 
This problem was already pointed out in Ref. JOj and is related to the well- 
known ordering problem in the quantization of classical Hamiltonians with position- 
dependent masses. Similar issues also arise in solid state physics Borrowing from 
the methods used to overcome these difficulties we use the following definition 
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where Iq = l a and In = h, and K(l, V ', y) is defined as: 
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In the continuum limit the partition function becomes 
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Figure 2. Phase diagrams for (a) rilling and (b) wetting transitions. The 
thick and dashed lines in both diagrams correspond to continuous and first-order 
boundaries between bound and unbound interfacial states, respectively. The 
arrows show representative paths along which continuous unbinding occur: (i) 
and (ii) for tricritical filling (critical wetting) and (iii) for critical filling (complete 
wetting), respectively. The filled circles represent the tricritical filling and critical 
wetting points, respectively. See text for explanation. 



where the complete orthonormal set satisfy 



( ad 


'1 d~ 


\ ~ 4S,dl 


Idl 



3a \ 



(7) 



In the thermodynamic limit Y — > oo we obtain the probability distribution function 
(PDF) for the midpoint interfacial height Pw{lo) — |^o('o)| 2 ) the wedge excess free 
energy fw — Eq and the longitudinal correlation length £ y = 1/(E\ — Eq), At this 
point we must remark that any definition of the path integral which is invariant upon 
exchanging l a and lb leads to an Schrodinger equation similar to Equation (JJJ) but with 
a different coefficient for the extra I /I 3 term in the effective binding potential |1U| . 

The change of variables A = ^/SZ/^l 3 / 2 /3 and t/>(0 = (2£Z/a) 1/4 0(A(O) El 
transforms the Equation J7J to: 

1 > ( V W [1{\)] - *) 0(A) = E<j>{\) (8) 



2 d\ 2 



72A 2 
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with 1(A) = ^3A/ y/8Ti/aj . In general, there will be an interfacial bound state at 
bulk coexistence for 8 — a if the strength of the small Iq attraction between the gas- 
liquid interface and the substrates, which we will denote as u, is greater than some 
value u c . Consequently, the filling transition is first-order if u > u c and critical if 
u < u c . Tricritical filling is observed when u — u c emerges as a new relevant field (in 
the renormalization-group sense) . If W(l) ~ ~a/l p + b/l q , different scenarios may arise 
as the range of the binding potential is varied. In particular, for p > 4 the long-range 
behaviour of Vw is dominated by the l/l 3 term for 9 = a, so the filling phenomena are 
fluctuation-dominated. This finding is consistent with the existence of two different 
fluctuation regimes for critical filling: mean-field if p < 4 and fluctuation-dominated 
regime if p > 4 In the fluctuation-dominated regime, examination of Equation 
(JSJ shows that there is an analogy between 3D continuous filling and 2D continuous 
wetting (see Figure EJ, where the role at wetting of the bulk ordering field h and the 
potential strength w are played by 9 — a and u for filling phenomena, respectively. In 
particular, 3D tricritical (critical) filling is analogous to 2D critical (complete) wetting, 
respectively. Different critical exponents which characterize the divergence of length 
scales can be defined. In addition to the critical filling critical exponents (3w and u y 
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defined along route (ii) in Figure^a) (see above), we can define new critical exponents 
for tricritical filling at 9 = a (route (i) in Figure |2fa)) as: 

i w ~ (u- u c y « , £, y ~(u~- u c y v l (9) 

and £j_ ~ £,y w , where £ w is the tricritical wandering exponent which in general may be 
different from C,w ( m contrast with the wetting case). More generally, in the vicinity of 
the tricritical point we anticipate scaling e.g. £ y ~ \u — u c | — "vA [(# — a)\u — u c | _A ] 
with the gap exponent A*. Thus along route (ii) £ v ~ (9 — a)" v yl^ . 

We focus now on the case of short-ranged forces as the prototype of the 
fluctuation-dominated regime. In addition to the hard- wall condition, Viy(^o) can 
be modelled as a contact-like attraction with strength u. An analysis of Equation 
for l — ► shows that the short-distance expansion of the PDF is either P\y ~ l n 
or Pyy ~ Zq. We anticipate that the former corresponds to tricritical behaviour 
and the latter to critical filling. It is remarkable that thermodynamic consistency 
at critical filling is ensured as the local density at the wedge bottom is non-singular, 
i.e. p w (0) — pi ~ T — Tf, where pi is the bulk liquid density 14 . This property is only 
obtained if the partition function is defined by Equation Q and Equation (JSJ). Thus, 
the ambiguity in its definition can be removed by imposing this regularity condition 
on the short-distance expansion of the interfacial height PDF. 

We report now our explicit results (details will be presented elsewhere). Along 
route (i) we find that there is only one bound solution to Equation Q for u > u c m 
1.358 with E oc (u — u c ) 3 and associated PDF 



Pw{k) 
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An, 

where Ai(x) is the Airy function and in the scaling limit ~ \u — m c | _1 - Thus 
lw ~ £± ( u — Uc)^ 1 and £ y oc (u — u c )~ 3 identifying (3 W = 1, v* = 3 and obtaining 
Cw ~ Civ = 1/3. As predicted, the short-distance behaviour of the PDF is linear with 
io. 

On the other hand, the scaling of the PDF for 9 > a is given by (see also Figure[3J): 



PwQo) <x 'oexp 



2e/o 2/q 



H, 
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where & = X" 1/2 [(6> / 'a) 2 - 1]~ 1/4 , e = £E $/a, v = e 2 /4 - 1/2 and H v (x) is the 
Hermite function The value of e is obtained as the smallest solution of the 

following equation: 



where the positive (negative) sign corresponds to u > u c (u < u c ), respectively. 
The inset of Figure plots the solution of this equation. As anticipated, note that 
scaling is obeyed in the vicinity of the tricritical point as the wedge excess free energy 
fw ~ F±{£e/£u)- For u> u c and £e/£ u >■ 1, we have checked numerically that the 
PDF H11J) converges to the expression given by Equation (|1C)|) for the corresponding 
value of £ u given by Equation (|12(l . Thus, the interface remains bound to the substrate 
when 9 — > a, in agreement to the first-order character of the filling transition. The 
thermodynamic path (ii) to the tricritical point corresponds to £ u — > oo, which 
corresponds to e w 1.086. Thus along this route lw ~ oc (9 — a) -1 / 4 similar 



Tricritical w 




Figure 3. Plot of the scaled PDF for e = —1.5 (thick dashed line), and along 
routes (ii) and (iii) in Figure a*), i.e. for e PS 1.086 (thick continuous line) and 
e ?3 1.639 (thick dot-dashed line), respectively. For comparison, the PDF from 
Equation I1UI with £ u as 1.968^9 (which corresponds to e = —1.5, see inset) is 
also plotted (thin dashed line). Finally, the scaled PDF obtained in Ref. 1101 
is also shown (thin dot-dashed line). Inset: Plot of e as a function of £,e/£u f° r 
u < u c (continuous line) and u > u c (dashed line). 

to critical filling. From analysis of the spectrum it is also possible to show that 
£ y oc (9 — a)~ 3 / 4 , so the tricritical gap exponent A* = 4. Finally, for thermodynamic 
paths (iii) far from the tricritical point, i.e. £ u — > 0, we have found that the scaling of 
the PDF is of the form shown in Equation l|llfl with e rj 1.639. Note that Pw ~ as 
/o — >• 0, in agreement with our previous statement. This condition is not fullfilled by 
the solution presented in Ref. JUj , although globally it does not differ too much from 
our exact solution (see Figure EJ). 

We finish by mentioning a remarkable connection for short-ranged forces between 
3D wedge filling and 2D wetting with random-bond disorder |15| . The critical 
exponents corresponding to tricritical and critical wedge filling can be obtained from 
generalized random- walk methods |16j in terms of the wedge wandering exponent £w 
In particular, they are found to have the same dependence of the critical exponents for 
critical and complete wetting, respectively, but in terms of an effective 2D wandering 
exponent equal to 2(w F° r short-ranged forces ((w = 1/3), this implies that the 
set of critical exponents is the same as for 2D random-bond disorder [TJ\. These 
predictions may certainly be tested in Ising model simulation studies and would be a 
stringent test of the theory of 3D wedge filling. 
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